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1 Introduction 

Film cooling injection is widely applied in the thermal design of turbomachinery, as it 
contributes to achieve higher operating temperature conditions of modern gas turbines, and 
to meet the requirements for reliability and life cycles. It is a significant part of the high- 
pressure turbine system. The film cooling injection, however, interacts with the main flow and 
is susceptible to have an influence on the aerodynamic performance of the cooled components, 
and through that may cause a penalty on the overall efficiency of the gas turbine. The main 
reasons are the loss of total pressure resulting from mixing the cooling air with the mainstream 
and the reduction of the gas stagnation temperature at the exit of the combustion chamber to 
a lower value at the exit of nozzle guide vane. In addition, the impact of the injected air on the 
evolution of the trace species of the hot gas is not yet quite clear. This work computationally 
investigates the film cooling influence on post-combustor trace chemistry, as trace species in 
aircraft exhaust affect climate and ozone. 

2 Codes used 

CGLENN-HT and NCC are selected as the primary working codes, while CNEWT-GRC is 
used as a reference. A description of these codes can be found in Ref. [1] 

3 Geometry and mesh 

The turbine vane of this study is based on an Allied-Signal film-cooled engine design [2]. 
An enlarged view of the cooling hole meshes in the leading edge region is shown in Figure 1. 
In the present study, a two-dimensional representation of the cooling holes is used, i.e. a slot 
is used to replace a row of holes. The plenums are also not included in the computational 
domain. 

An overset /Chimera grid is generated for both CGLENN-HT and NCC. This grid consists 
of two component grids: a hyperbolic marching grid generator in T3D is used to generate an 
O grid around the vane, while an algebraic H grid is generated to cover the flow passage. 
OVERGC[3], acronym of overset grid communicator, is then used to build grid interfaces 
between these two component grids. Furthermore, FPLOT is used to convert the composite 
grids resulting from OVERGC to the format of overset unstructured grids for NCC, as well as 
the format of overset structured grids for CGLENN-HT. FPLOT is also used to tag boundary 
elements and fluid elements for the NCC. The final composite grid is shown in Figure 2. The 
positions and directions of the cooling air slots are indicated in Figure 3 
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Figure 1: The leading edge region of a three dimensional mesh of the cooling holes 


A TURBINE VAIME WITH 12 COOLING SLOTS 



Figure 2: A two-block overset grid around the vane. (Thickness of vane is about .01 meter.) 
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Figure 3: The positions and directions of the cooling air slots. 


Slot no. 

M ( kg \ 
A v s-m 2 ' 

flow angle (degree) 

i 

110.46 

-65.5 

2 

97.38 

-62.9 

3 

130.54 

-63.8 

4 

115.52 

-61.4 

5 

94.52 

186.2 

6 

186.28 

182.7 

7 

202.12 

180.9 

8 

210.6 

179.1 

9 

251.68 

177.1 

10 

128.76 

175.6 

11 

159.66 

44.5 

12 

171.26 

47.0 


Table 1: Film-cooling inlet flow parameters 
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4 Boundary conditions 


At the mainstream inlet, the total pressure of the mainstream hot gas is set to be 1629080 
N/M 2 ; the total temperature is 1644.42 K ; and the specific heat at constant pressure is 
1138 J/kgK, Following standard practice the value for the turbulence intensity is 4% of the 
computed inlet normal velocity, while the turbulence length scale at the inlet is set to be 
.05 m . At the exit, the static pressure is set to be 938350.08 N/M 2 , i.e. 57.6% of the inlet 
total pressure. The surface temperature of the vane is 1151 K . 

The coolant mass flow is about five percent of the mainstream hot gas mass flow. Cooling air 
at a total temperature of 822.1 K is discharged into the mainstream through 12 slots along the 
vane surface to form a cooling film. The specifications of the coolant mass flux and discharged 
angle at each slot are listed in Table 1. The composition of the species for the inlet mainstream 
gas as well as the cooling air are listed in Table 2. The trace chemistry mechanism used is the 
one reported in [4]. It has twenty five species and seventy four reation steps. 



Species 

Mass fraction at main inlet 

Mass fraction of cooling air 

1 

C(S) 

.2437335E-05 

0 

2 

H2 

.2327455E-07 

0 

3 

02 

.1438427E+00 

0.21 

4 

H20 

.2986734E-01 

0.01 

5 

O 

.1098798E-05 

0 

6 

H 

.4135477E-09 

0 

7 

OH 

.4075543E-04 

0 

8 

H02 

.7921337E-06 

0 

9 

H202 

.2097706E-07 

0 

10 

CO 

.1947870E-03 

0 

11 

C02 

.7671610E-01 

0.02 

12 

N 

•0000000E+00 

0 

13 

N2 

.7491807E+00 

0.76 

14 

NO 

.1098798E-03 

0 

15 

N02 

.1877946E-04 

0 

16 

HNO 

.1048853E-08 

0 

17 

HONO 

.1388482E-06 

0 

18 

N20 

•0000000E+00 

0 

19 

N03 

.474481 IE-10 

0 

20 

HN03 

.4505073E-09 

0 

21 

SO 

.3026690E-10 

0 

22 

S02 

.2377400E-04 

0 

23 

S03 

.6762604E-06 

0 

24 

HS03 

.1678165E-09 

0 

25 

H2S04 

.2727018E-10 

0 


Table 2: Species distributions at the inlet of the and cooling air 
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5 Evaluation of transport properties 


Transport properties are evaluated differently in different codes, they are summarized as 
the following three options. 

1. Option 1 — The temperature is iteratively computed from the static enthalpy of the gas 
mixture, /i, and the specific heat at constant pressure of the gas mixture, C p . In addition, 
the gas constant of the mixture, the specific heat at constant pressure of the gas mixture, 
laminar viscosity of the mixture and thermal conductivity of the mixture are all functions 
of species mass fractions and the temperature. 

2. Option 2 — The temperature is computed from the enthalpy of the gas mixture, h and 
a constant C p . Gas constant of the mixture, specific heat at constant pressure of the gas 
mixture, C pi laminar viscosity of the mixture and thermal conductivity of the mixture are 
all constants which are either user specified parameters or computed via other parameters 
such as Prandtl number. 

3. Option 3 — The gas constant of the mixture and the specific heat at constant pressure 
of the gas mixture, C pi are constants provided by the users. The laminar viscosity of 
the mixture is calculated using a power- law for its dependence on temperature.. The 
thermal conductivity of the mixture is computed through Prandtl number and C p . The 
temperature is computed from the equation of state. 


6 Computational results 


Five sets of computed results are assembled together to show the difference due to various 
codes and cooling air injection. They are: (1) CGLENN-HT using option 3 without cooling 
air, (2) NCC using option 1 without cooling air, (3) CNEWT-GRC using option 2 without 
cooling air, (4) CGLENN-HT using option 3 with cooling air, and (5) NCC using option 1 with 
cooling air. Numerical simulation without film cooling represents the first step in the analysis. 
Its converged solution is then used as the initial condition for the film cooling computation. 
Data reduction has been performed to get one dimensional distribution along the axis of the 
flow path by averaging the CFD results. 

The pressure contours of the five cases are shown in Figure 4. The area-averaged pressure 
distributions along the axis of the flow path are given in Figure 5. These results suggest that 
the differences of the pressure distribution among the five cases are minimal. 
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CNEWT 

Option 2 

WITHOUT COOLING AIR 


NCC 

Option 1 

WITHOUT COOLING AIR 


C GLENN -HI 

Option 3 

WITHOUT COOLING AIR 



PRESSURE 
11 1.6415Se+0B 
10 1 , 539230+06 

H I ,430S7e+Q6 
1 33451e+06 
7 1 .2321 50+06 
$ 1.12979S+06 
5 1 0£743e+06 
4 925071 
3 622711 
| 720351 
| 017991 



NCC 
Option 1 
COOLING AIR = 


CGLENN-HT 
Option 3 

.05 INL ET AIR COOLING AIR = 


.05 INLET AIR 


i 


Figure 4: Comparison of pressure distributions. 


Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 

Area- Averaged Pressure Distribution Along Axis of Flow Path 



Figure 5: Comparison of area- averaged pressure distributions along the axis of flow path. 
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CNEWT 

Option 2 

WITHOUT COOLING AIR 


NCC 

Option 1 

WITHOUT COOLING AIR 


CGLENN-HT 

Option 3 

WITHOUT COOLING AIR 



Figure 6: Comparison of temperature distributions. 


The temperature contours of the results are presented in Figure 6. The temperature dis- 
tributions along the axis of flow path are depicted in Figure 7. It indicates that the impact 
of film cooling on the temperature is significant. It also suggests that the evaluation of the 
transport properties affects the temperature more for the case with cooling air than the case 
without cooling air. 


Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 


Area- Averaged Temperature Distribution Along Axis of Flow Path 



Figure 7: Comparison of area- averaged temperature distributions along the axis of flow path. 
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CO distributions are presented in Figure 8 and Figure 9. The mass fraction decreases as 
the flow passes the turbine vane with or without cooling air. 


CNEWT 

Option 2 

WITHOUT COOLING AIR 


NCC 

Option 1 

WITHOUT COOLING AIR 


C GLENN -HI 

Option 3 

WITHOUT COOLING AIR 



MASS FRACTION 


Figure 8: Comparison of CO distributions. 


Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 


Area- Averaged CO Distribution Along Axis of Flow Path 



Figure 9: Comparison of area-averaged CO distributions along the axis of flow path. 
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NO distributions are shown in Figure 10 and Figure 11. The mass fraction for the case 
with cooling air decreases significantly as the flow passes the turbine vane. 


CNEWT 

Option 2 

WITHOUT COOLING AIR 


NCC 

Option 1 

WITHOUT COOLING AIR 


C GLENN -HI 

Option 3 

WITHOUT COOLING AIR 



NO 

11 0.0001 24&S9 
m 3.0001 12103 

M 9.9B472e-Q5 
S,71915e-05 
f 7.47354Q-05 
& 6.227&5e-05 
5 4&02^a-O5 
4 3.73677S-05 
3 2 491 13*- 05 
2 1 .24559e-Q5 



MASS FRACTION 



J 

NCC 

Option 1 
COOLING AIR = 


CGLENN-HT 

Option 3 

,05 INLET AIR COOLING AIR - 


.05 INLET AIR 


Figure 10: Comparison of NO distributions. 


Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 

Area- Averaged NO Distribution Along Axis of Flow Path 



Figure 11: Comparison of area- averaged NO distributions along the axis of flow path. 
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N0 2 distributions are given in Figure 12 and Figure 13. The mass fraction with cooling 
air increases more than that without cooling air as the flow passes the turbine vane. 


CNEWT 

Option 2 

WITHOUT COOLING AIR 


NCC 

Option 1 

WITHOUT COOLING AIR 


C GLENN -HI 

Option 3 

WITHOUT COOLING AIR 




N02 

11 2.8931 5e-Q5 
Id 2 .603330-05 

■ 2.314620-05 
|z.Ci252e-05i 
7 1 .735390-05 
6 1 446579-05 
5 1 , 157269-05 
4 0.67345#-M 
3 5.70639-06 
2 2.893150-06 

MASS FRACTION 


CGLENN-HT 

Option 3 

INLET AIR COOLING AIR - .05 INLET AIR 


Figure 12: Comparison of N0 2 distributions. 


Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 

Area- Averaged N02 Distribution Along Axis of Flow Path 



Figure 13: Comparison of area-averaged N0 2 distributions along the axis of flow path. 
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SO 2 distributions are shown in Figure 14 and Figure 15. The mass fraction with cooling 
air decreases by an amount of 25 to 30 percent of the inlet value as the flow passes the turbine 
vane, while the mass fraction without cooling air remains essentially the same. 


CNEWT 

Option 2 

WITHOUT COOLING AIR 


NCC 

Option 1 

WITHOUT COOLING AIR 


CGLENN-HT 

Option 3 

WITHOUT COOLING AIR 


S02 

11 2.38195* -05 
IQ 2.14376e-05 

I I ,90556e-05 
1 ,66737e-05 

1.429176-05 

6 1.19098e-05 
5 9, 52786-08 
4 7.145856-06 
3 4.76396-08 
2 2.381 95e-Q6 

MASS FRACTION 



Figure 14: Comparison of SO 2 distributions. 


Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 


Area-Averaged S02 Distribution Along Axis of Flow Path 



Figure 15: Comparison of area-averaged SO 2 distribution along the axis of flow path. 
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SO% distributions are presented in Figure 16 and Figure 17. The mass fraction with cooling 
air increases more than that without cooling air as the flow passes the turbine vane. 


CNEWT 

Option 2 

WITHOUT COOLING AIR 


NCC 

Option 1 

WITHOUT COOLING AIR 


C GLENN -HI 

Option 3 

WITHOUT COOLING AIR 



£03 

11 1.70925e-O6 
10 1 .538336-0 6 

M l ,3G74fl-06 
1.19646Q-06 
7 1 02555a-06 
& 0 

5 S 8S7e-07 
4 5.12775*-07 
3 3.413&MD7 
2 1 .709250-07 



MASS FRACTION 


Figure 16: Comparison of SO% distributions. 



4 

NCC 

Option 1 

COOLING AIR - 


CGLENN-HT 

Option 3 

.05 INLET AIR COOLING AIR = 


.05 INLET AIR 


Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 


Area- Averaged S03 Distribution Along Axis of Flow Path 



Figure 17: Comparison of area-averaged SO% distribution along the axis of flow path. 
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The area- averaged HONO distributions are shown in Figure 18. The mass fraction with 
cooling air varies more than that without cooling air. 


Effect of Film Cooling of A Generic Turbine Vane 


Area- Averaged HONO Distribution Along Axis of Flow Path 



Figure 18: Comparison of area-averaged HONO distributions along the axis of flow path. 
Figure 19 gives the area-averaged velocity distribution. 

Trace Chemistry Effect on Film Cooling of A Generic Turbine Vane 


Area- Averaged Velocity Distribution Along Axis of Flow Path 



Figure 19: Comparison of area-averaged velocity distributions along the axis of flow path. 
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7 Concluding Remarks 

• A two-dimensional representation of the 3D geometry is used for the present analysis. 

• The coolant mass flow rate is about five percent of the inlet main flow. 

• The present work indicates that the cooling air has an influence on the temperature and 
trace species mass fractions. In particular, the values of the mass fractions of NO and 
SO 2 are lower while the mass fractions of NO 2 and SO 3 are higher than the correspoding 
values without the presence of the cooling air. 

• Among other factors such as numerical schemes and turbulence models, a significant 
factor leading to the difference of the results obtained by NCC and CGLENN-HT is the 
approach for evaluating the transport properties. The effects of the heat release due to 
the chemical reactions are not considered in the CGLENN-HT code. 

• The CGLENN-HT code was designed to run in parallel on a cluster of 2, 3, 4, 5, 6 or 
12 nodes using MPICH, a freely available implementation of the MPI (message passing 
interface) standard. Only the chemistry subroutines of CGLENN-HT are parallelized. 
The NCC code was designed to run in parallel on a cluster of arbitray number of nodes 
using MPICH. The entire code is parallelized based upon the principle of the domain 
decomposition. Thus wall clock time of computation using the NCC code can be sig- 
nificandly reduced by requesting more computer nodes, while the CGLENN-HT is much 
more limited in this respect. For the film cooling cases running NCC and CLGENN-HT 
on the LOMAX machine of the NAS facility, the run time statistics is 

— NCC — On a single node, the average wall clock time is about 120 seccond per 
iteration. It can achieve a speedup in an almost linear fashion by adding more 
CPU’s. For example, on 32 nodes, the average wall clock time is about 3.75 second 
per iteration. 

— CLGENN-HT — On a single node, the average wall clock time is about 80 second 
per iteration. The average wall clock time is about 59 second per iteration on two 
nodes, 46 second per iteration on three nodes, 41 second per iteration on four nodes 
and 39 second per iteration on six nodes. 

This speedup limitation of the CGLENN-HT code is taken into consideration, as we decide 
to use NCC for furthering the post-combustion chemistry modeling and simulations. 
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